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In any realistic particle physics model of inflation, the inflaton can be expected to couple to 
other fields. We consider a model with a dilaton-like coupling between a U(l) gauge field and 
a scalar infiaton. We show that this couphng can result in observable non-gaussianity, even in 
the conventional regime where inflation is supported by a single scalar slowly rolling on a smooth 
potential: the time dependent inflaton condensate leads to amplification of the large-scale gauge field 
fiuctuations, which can feed-back into the scalar/tensor cosmological perturbations. In the squeezed 
limit, the resulting bispectrum is close to the local one, but it shows a sizable and characteristic 
quadrupolar dependence on the angle between the shorter and the larger modes in the correlation. 
Observable non-gaussianity is obtained in a regime where perturbation theory is under control. If 
the gauge field is identified with the electromagnetic field, the model that we study is a realization 
of the magnetogenesis idea originally proposed by Ratra, and widely studied. This identification 
(which is not necessary for the non-gaussianity production) is however problematic in light of a 
strong coupling problem already noted in the literature. 
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I. INTRODUCTION 



Over the next few years Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) probes will 
measure the primordial cosmological fluctuations with substantial improvements in accuracy and over a widening 
range of scales. With this wealth of data there is an exciting prospect to strongly constrain, or perhaps measure, 
nongaussian statistics of the primordial density perturbations. Nongaussian effects encode a wealth of information 
about the physics of the very early universe and might provide a powerful tool to discriminate between different 
models. (See for a recent review.) 

It is often claimed that primordial nongaussianity will be undetectable in the simplest models, where the inflation is 
driven by a single field that is slowly rolling along a smooth, flat potential. The physical reason is that nongaussianity 
is a measure of the strength of interactions, while the requirement of a flat potential usually constrains inflaton self- 
interactions to be weak ■ A number of models have been constructed which do produce a detectable noiigaussian 
signature, for example using sound speed effects Q, higher derivatives 0], non- vacuum initial conditions 
sharp potential features [11, jL2], post-inflationary effects [13, [3i ^tc. 

However, in general there are not only inflaton self-couplings, or gravitationally suppressed couplings. For instance, 
in any realistic particle physics framework, the inflaton field ip can be expected to couple to "matter" fields. "'^ Such 
couplings are certainly necessary for successful reheating after inflation, and they are much less constrained by the 
requirement of slow roll. Their consistent inclusion can radically impact the phenomenology of the model; see the 
works p"5l - l23 | for explicit examples and scenarios. 

We focus our attention to inflaton interaction with gauge fields. For a singlet inflaton, there are two very natural 
classes of gauge field interactions to consider, depending on the parity. For a pseudoscalar inflaton one expects axial 
interactions of the type 



^pseudo _ _ ^ p^jy p /-.N 



where F^^ — df^Ai, — d,yA^ is the field strength associated to some U{1) gauge field and F'^'^ — e^""^ Fa/3/2 is its 
dual. For a scalar inflaton, on the other hand, one may expect couplings of the form 

/'Scalar 

•'-int ~ ^ ^ ^fJ-f , l^y 

where I{(p) plays the role of a field dependent gauge coupling. The interaction ^ is typical of moduli or dilaton-like 
fields in string theory and supergravity frameworks. As noted in psj . coupling the inflaton to a gauge field is actually 
the only way to reheat for some of these models. 

Pseudoscalar couplings of the type ([T]) have been studied extensively in [l7l - [l9j . For natural values of /, such 
couplings lead to a copious production of gauge field fiuctuations that feed-back on the scalar and tensor cosmological 
perturbations. Observable nongaussianity of the equilateral type is naturally generated, without any exotic model- 
building ingredients Moreover, this is correlated also with a gravitational wave signal that can be detectable 
with interferometers ^19, . .23i | . 

In this paper, we show that a similarly rich phenomenology is possible also in simple models with a single scalar 
inflaton in slow roll, via the coupling The underlying mechanism is quite novel. The time-dependence of the 
inflaton condensate ipo{t) = {'~p{t, x)) breaks the conformal invariance of the gauge field sector and leads to amplification 
of the quantum fluctuations of A^, similarly to the well-known mechanism that produces scale invariant curvature 
fluctuations during inflation. For simplicity, we focus our attention on the case where large scale fluctuations of the 
gauge field are produced during inflation with a scale-invariant "magnetic" component. We notice that the same 
coupling that leads to production of gauge field fiuctuations also implies that these produced fiuctuations must, in 
turn, couple to the cosmological perturbations of the infiaton, Sip(t,x) — ip{t,x) — tpoit)- We find that the feed-back 
of produced gauge fiuctuations on Sip contributes a new component to the observable curvature fiuctuations that is 
highly nongaussian and is uncorrelated with the usual spectrum from quantum vacuum fluctuations. For reasonable 
parameters, we obtain nearly local-type nongaussianity with shape 

1 -I- cos^ ffci, A-aj 

(Cfei Cfe, Cfe3 ) oc Tap + permutations (3) 

"-1 "-2 



Here we use the word "matter" to describe any fields that do not play a significant role in driving the infiationary expansion. 
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and at the level /nl = 0(10 — 100). Such values are close to current observational limits and will be probed in 
the near future. The shape ([3]) has a strong overlapping with the local template in the squeezed limit, where both 
shapes are enhanced. However, in this limit, (|3]) has a quadrupolar dependence on the angle between the shorter side 
and either of the longer sides. Contrary to what typically happens for non-gaussianity sourced by scalar fields [2^, 
this angular dependence does not vanish in the squeezed limit, and it contributes to about 1/5 of the amplitude of 
This therefore appears as a signature of non-gaussianity from higher spin fields, and it may be an important 
distinguishing feature when the model is confronted with observations. ^ 

Our scenario is a realization of the "feeder" mechanism [27\ : consequently, the Probability Distribution Function 
(PDF) exhibits a non-hierarchical structure. This unique feature could help to observationall y d istinguish the model 
from other constructions that give nearly local nongaussianity (for example the curvaton (28l, l29l|). 

Models of vector fields producing non-gaussianity were also proposed in [sOl . In these models the vector is amplified 
either from the time dependence of its kinetic term, as we consider here, or from a nonminimal coupling (as an effective 
mass term) to the scalar curvature R. This second case, however, introduces a longitudinal vector component, which 
turns out to be a ghost [HI, HI]. (In some of these models, the energy in the vector field sources anisotropic inflation; 
we do not study this possibility here, but we refer the reader to Rcf. [33] for a review.) In these models, the vector 
acts as a curvaton, and it contributes to the spectrum and the bispectrum through its energy density. We compute 
instead the non-gaussianity resulting from the same coupling that leads to production of gauge field fiuctuations, 
and which remains imprinted in the inflaton perturbations even when the energy density in the gauge field is very 
negligible at reheating. The model of ref. [M] studied non-gaussianity from a g auge field amplified analogously to 
what we do here, but imprinted through the waterfall field of hybrid inflaton [35j. 

It is interesting to contrast the mechanism considered in this paper to the closely related physics of inverse decay 
that was studied in [l3, [13 • In the case at hand, the feed-back of the produced gauge field fluctuations on the 
inflaton fiuctuations leads to a super-horizon growth of the curvature perturbation. Such growth is consistent, since 
the produced gauge field fiuctuations provide a source of large scale iso-curvature perturbations. Consequently, we 
find a bispectrum which is very close to the local shape. On the other hand, for axial couplings one finds [13] that the 
relevant production of inflaton perturbations arises near horizon crossing, and then the bispectrum is instead nearly 
equilateral. 

Models of the type which we study have received considerable attention in the literature, in connection with 
primordial magnetogenesis [36l - l45j . At face value, our choice of P{(p) can produce large scale magnetic flelds with 
a sufficient amplitude to account for observations at galaxy and cluster scales; see [46-50] for reviews. This would 
open the interesting possibility of correlating the magnetic field with the primordial perturbations [1^ , although the 
correlation would only involve the component of metric perturbations that are sourced by the vector field, and that 
is typically subdominant with respect to the vacuum part. Moreover, the magnetic field would induce non-gaussianiy 
from its direct coupling to the CMB photons [Sll. This effect can be observed for a magnetic field at the ^ nC— 10 nC 
level, while the non-gaussianity we have obtained arises from the direct coupling to the inflaton that generated the 
gauge field, and can be observable even if the current "magnetic" filed is significantly smaller. (A list of works that 
study the general signatures of a magnetic field on the CMB can be found in the review |48|.) 

One major problem with identifying the gauge field with the electromagnetic one is due to the fact that a scale 
invariant magnetic field can only be obtained if the effective gauge coupling, g{t) — I [ipo{t)]^^ , decreases by many 
orders of magnitude during infiation (specifically, by a factor of e^^'°' , where A'tot is the number of e- folds of infiation) . 
If one starts with from — 0(1), then the gauge coupling at the end of infiation will be much too small to identify 
Afj^ with the Standard Model (SM) photon. Normalizing instead the coupling constant to be the electromagnetic one 
at the end of inflation would instead entail an unacceptable breakdown of perturbation theory. This problem was 



stressed in [4j] , and it is a serious obstacle in identifying the gauge fleld of the mechanism with the photon. In the 



following, we discuss some unsuccessful attempts of solving this problem. We cannot of course rule out that a solution 
of the problem can be found, but we believe that it would require a substantial modification of the model. 

As a consistency check on our calculation, we have verified that perturbation theory is well under control using a 
variety of different diagnostics, including backreaction effects, the amplitude of curvature perturbations, the energy 
density in fiuctuations, and the Weyl tensor. We show that the "new" metric perturbations which are sourced by the 
gauge field fluctuations are generically sub-dominant, as compared to the standard vacuum contribution. 

The organization of this paper is as follows. In Section we introduce our model and compute the production 
of gauge field fluctuations. We then discuss the challenging connection to magnetogenesis. In Section IIIII we give a 
succinct review of the key phenomenological predictions for the spectrum and bispectrum of curvature fluctuations, 
and also the spectrum of tensor fiuctuations. In Section IIVI we present a detailed computation of the 2-point and 



^ We thank Antony Lewis for stressing the importance of this in a private communication. 
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3-point correlation functions of the scalar perturbations. In Section |V] we present the computation of the 2-point 
correlation function of gravitational wave fluctuations. In Section I VII we discuss the validity of our perturbative 
analysis. Finally, in Section IVIIi we conclude. Appendix A gives some technical details on second order cosmological 
perturbation theory in the spatially flat gauge. In this paper we compute the feed-back of the produced gauge field 
fluctuations on the inflaton using the Green function method that was developed in ^20l . [22j and employed also in 
This formalism is equivalent to the "in- in" approach at leading order. In Appendix B we demonstrate this 
equivalence explicitly for the case at hand. (See also |19|.) 



II. GAUGE FIELD PRODUCTION 



The Model 



We consider a simple model with a dilaton-like coupling of the inflaton to a U(l) gauge field 



S = d X\/—g 



Mi 1 , 



(4) 



In this action, F^^ is the field strength of A^, and / and V are functions of the inflaton (p. We assume that the 
potential V{(p) is sufficiently flat to support a long phase of quasi de Sitter expansion. As usual, we require that 
e, I77I <^ 1 where the slow roll parameters are 



/V. 



2 \ V ' P V 



(5) 



To simplify the calculation, we choose I{(p) and V{ip) to be related to each other in such a way as to obtain the 
solution I (X oT" , where a is the scale factor of the Universe, during inflation. The required relation can be obtained 
[53I I by taking the ratio of the slow roll relations 



1 



3M2 



V , 3 Hip 



(6) 



(assuming that the interaction with the gauge field provides negligible corrections to the dynamics, see below). Here 
and in the following, dot denotes derivative with respect to physical time, ", y?" denotes derivative with respect to 
ip, while H is the Hubble rate, H — ^. In this way, one forms a differential equation for da/ dip, that is integrated 

The functional form of / (ip) can be then set to the n — th power of the right hand side 



into a oc exp 



expression. For definiteness, we consider a monomial inflaton potential 



V ^ '■93'' , / = /end exp 



nip 
2rM2 



(7) 



although other choices are clearly possible. Here /end is the value of the coupling function at the end of inflation, 
when the inflaton is in the vacuum ip = Q. We therefore assume that, after inflation, the function / sets to a constant. 
We could then normalize /ond = 1- 

We found it algebraically convenient to deflne "electric" and "magnetic" components of the gauge field as 



B - r) A 



(8) 



where here and in the remainder of this work the Coulomb gauge = is assumed. We do not necessarily assume 
Af^ is the Standard Model photon (more on this later), however, we will sometimes use the language "electric field" 
and "magnetic field" , by analogy with standard electromagnetism. 

As the gauge field has no classical expectation value, its perturbations do not couple to that of the inflaton or of 
the geometry at linearized order. We can therefore solve for these perturbations by assuming a FRW background, 
and by treating / as a classical function. Therefore, we can simply set / cx a" for the remainder of this Section. The 
time component of the vector equation of motion is solved by diAi — 0, and we can decompose the vector potential 
as 



A 



E 

A=± 



d^k 



(27r) 



3/2 



e-A (fc) e*'"-^^ [oA (k) Ax (k) + a{ (-fc) Al (fc) 



(9) 
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where the circular polarization operators satisfy k- £± (^k^ — 0, kx e± (jtj = ^ike± (jt^ , e± (~^') ^ ^± (^) ' ^'^"^ 

normalized according to ex ^fc^ • ex' (jtj = Sxx'- The annihilation / creation operators satisfy ax (jtj , a^, ^A/^ = 

6xx' 5(3) {k-k'). 

The vector mode functions then satisfy 

Vi' + (^k^ - Vx ^ , Vx=IAx (10) 

where prime denotes derivative with respect to confornial time r. For a constant /, one recovers the typical Minkowski 
e~''^'^/v^ solution for the mode functions, due to the fact that the gauge field is conformally coupled to the FRW 
metric. For our purposes, it is sufficient to obtain the leading expression of the mode functions in a slow roll expansion. 
Namely, for the de Sitter limit a — —7^, the mode solution that, up to an arbitrary constant phase, is normalized to 
e~*'^'^/\/2fc (the so called adiabatic vacuum) in the asymptotic early time / high momentum regime is 

Vx = ^y^V^Hl%^{-kr) (11) 

This solution has been discussed at length in (43 | for all values of n, and it is unnecessary to review all their study 
here. We only discuss the n = 2 case, which results in a scale invariant "magnetic" field. Interesting non-gaussian 
properties of the primordial perturbations may be possible also for other values of n. The arbitrary phase in (ITT|) has 
been chosen so that the function Vx is real and positive in the super horizon limit. 

In this case, the "electric" and "magnetic" gauge field operators during inflation reduce to 



E-Ij07^e^'-'Ej: , Ej: = -i^j:,^ex{k)[ax{k)+a 
B-Ij07. , Bj:=^-^j:,X^ex (k) [ax (fc) + a] 



-kr < 1 (12) 



in the super-horizon limit. As for the standard scalar case, the mode function does not oscillate in the super-horizon 
regime, which is a signal that the field has become classical [s^l. The energy densities are given by 

PE='-^^^ jdkk{l + k-r-) , PB-^^^Jj, (13) 

We must only compute the energy in the classical fields, and we therefore limit the integrals to momenta that have 
exited the horizon during inflation. Modes of smaller wavelength remain in their vacuum state and their contribution 
to the energy density must be renormalized away (this is part of the cosmological constant problem). The smallest 
momentum fcmin — 37— corresponds to modes that have exited the horizon at the start of inflation; for any moment r, 
the largest momentum fcmax — 37 corresponds to modes that have just exited the horizon at that moment. Therefore, 
for any r ^ Tin during inflation, 

IbTT^ 47r^ a (rin) 

These behaviors are very different from the usual behavior p oc a^** for radiation. This shows that energy is being 
transferred from the inflaton to the gauge field through the I^F^ coupling. We also note that the energy in the 
"magnetic" component is greater than that in the "electric" one. The energy in the "magnetic" component is scale 
invariant, and the logarithmic increase in the final result is due to the increase of the phase space of the modes that 
have become classical. We also note, that, despite for most of the super-horizon modes the density is several orders 
of magnitude smaller than and decreases with time, the total value of pE is "only" logarithmically suppressed 
with respect to that of ps- This is due to the fact that the integral for pE has most of its support in the UV region, 
where the "electric" and "magnetic" energy densities are not too different from each other. 

In passing, notice that our choice to produce scale-invariant "magnetic" fields - as opposed to "electric" fields ~ is 
essentially arbitrary from the perspective of primordial nongaussianity. Indeed, there is an electric/magnetic duality 
that leaves the Maxwell equations invariant under the replacement E —B, B ^ E and / — >■ 1/J; see [s^l for more 
discussion. In this case at hand, this means that we can interchange the "electric" and "magnetic" spectra simply by 
taking n — ~2 rather than n = 2. The feed-back of these produced fluctuations on the scalar inflaton is essentially 
unchanged under such a replacement. 
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B. Backreaction Bounds 



Throughout the discussion above, we have assumed that the produced gauge field fluctuations have a negUgible 
effect on the homogeneous baclcground dynamics during inflation. To ensure that this assumption is consistent, we 
must verify several backreaction constraints. First, we require that the energy density in the "magnetic" field is much 
smaller than the potential energy driving infiation. Hence, we require 

TT — "T^ T7^ In — < 1 (15) 

V 47r2 M2 a (Tin) 

A second constraint arises because the produced gauge field fiuctuations modify the homogeneous Klein-Gordon 
equation for the infiaton condensate (fioit). To ensure that the usual slow roll equations ([6]) are reliable, we require 
that the "driving force" in the inflaton equation of motion is dominated by the derivative of the inflaton potential. 
That is, we require: 



\V^\ » 



(16) 



Finally, we note that the backreaction of produced gauge flelds can also lead to a correction for the effective 
mass of the inflaton. The easiest way to see this effect is to note that the action ^ contains a term of the form 
{l'^) ,vv{^'-pY {F'^) once we expand t/s = tpo + 5</5 and replace which its vacuum average, to estimate the magnitude 
of backreaction effects. We require that this "new" correction to the inflaton mass is much smaller than the Hubble 
scale so that we do not spoil the scale invariance of the spectrum. This amounts to a constraint 

^ 'B^) < (17) 



eM2 



The conditions (|16p and ([T7|) are more stringent than dTSl) . In both cases, we obtain a constraint on the total 
number of e-foldings of inflation of the form 

iVtot < IQ-^V-^ (18) 

where iVtot = In [ocnd/ain] and we have defined 

ttA tt2 

P ~ ~ -^—^ ~ 2.5 • 10-9 (19) 

which gives the amplitude of the power spectrum from the usual vacuum fluctuations ^55"] . (In this model there is also 
an additional contribution to the power spectrum coming from second order effects, however, we will always work in 
a regime where this is subdominant. More on this later.) 

The condition ((T5)) then indicates that backreaction is negligible provided that iVtot <C O (10''). Note that iVtot 
enters in this condition because the "magnetic" fleld energy grows (proportionally to the number of e-folds) during 
the whole duration of inflation. We stress that A'tot may be much greater than the number of e-folds Ncmb — 50 — 60 
which separates the moment at which the largest CMB scales exited the horizon to the end of inflation. 



C. Connection with Magnetogenesis 

Magnetic fields have been observed at many scales. They are present in structures (eg - galaxies, galaxy clusters and 
high redshift protogalactic structures) with strength ~ 10~^ — 10^^ G, and in the low density intergalactic space with 
strength ^ 10-^** — 10"^^ G. See j46i - t48i . i50i] for reviews. The origin of these fields is not well understood. Although 
an astrophysical mechanism is not ruled out, the observed homogeneity and large coherence length (~' kpc — Mpc) 
could suggest a primordial origin. 

For a standard electromagnetic action, the photon is conformally coupled to a FRW geometry; loosely speaking, 
the scale factor drops from the action term y^—gF^, and the photon remains in its vacuum state. Mechanisms 
for generation of magnetic field during infiation break the conformal invariance by introducing some extra-term. For 
instance, in [H^ couplings to the curvature invariants of the type RA^ and R^j^i^A'^A'^ were considered. These couplings 
however break the U(l) invariance associated to the electromagnetic field, and one should worry about the longitudinal 
photon component that they introduce. It was shown in [Slllsll that, with the R^A^ coupling introduced in [5^, the 
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longitudinal photon is a ghost. It is safer to consider models that preserve the U(l) invariance. Axial couplings j(pFF 
have been considered in [57l - l59| . In such models it is typically difhcult to produce a sufficiently large field. Note that 
any attempt to generate primordial magnetic fields via an axial coupling must take into account the limits on / due 
to nongaussianity from inverse decay effects [l7l - [l9l |. which are much more stringent than backreaction bounds. 

The model (|4|) has been studied in connection with primordial magnetogenesis (sgI - IssI . liol . [43l - l45l | . It is indeed 
tempting to identify with the standard model photon. If we do so, and we assume that the electromagnetic energy 
density scales in the standard way p cx from the end of inflation on, we find 

where we have assumed matter domination due to the coherent infiaton oscillations until reheating takes place at 
the temperature T^h- We have also disregarded the current departure from matter domination (this gives a negli- 
gible correction to the estimate), and treated the value of H as constant during inflation. We note that a lower 



reheating temperature results in a smaller value of the magnetic fleld today [44| : the estimate obtained in [45| as- 
sumes that radiation domination starts immediately after inflaton. In this case the expression (j20p evaluates to 
10-10 G(H/10i5GeV). 

The problem arising in associating the fleld with the electromagnetic photon has already been stressed in (4^. 
The model (jj]) must be supplemented by the action for the matter flelds. The most minimal approach is to assume 
that P {(f) only enters in the term, so that the relevant terms for the electromagnetic coupling of the (Standard 
Model) fermions are 

Cr..ttcr = -^-^F^ -i,j^'{d^ + ieA^),p (21) 
If this is the case, the "instantaneous" electric coupling constant is 

ephysical = e/-l [(y9o(t)] (22) 

During inflation we have / oc a" and, consequently, for n > 0, the electric coupling constant decreases by a huge 
factor during inflation (we recall that the scale invariant B fleld is obtained for n = 2). Thus, if we take ein ^ C>{1) 
at the start of inflation, then the gauge coupling at the end of inflation will be extremely tiny. Assuming no further 
evolution of Cphysicai (t) in the post-inflationary epoch, we clearly cannot identify with our photon. Alternatively, 
if we normalize / such that Cphysicai after inflation coincides with the present value, we necessarily imply that ephysicai 
was extremely large all throughout inflation, apart from the very last stages. Even if there were no real charged 
particles during inflation, this would lead to strong quantum effects from the vacuum fluctuations of these flelds, and 
to a quantum theory (at the very least) out of computational control. This poses serious questions on any result 
obtained from the model. We stress that this problem is not present if is a hidden sector gauge field, since in this 
case one may assume that its associated physical coupling constant is <, 0(1) at the start of inflation. 

We briefly comment on a few (unsuccessful) attempts to solve this problem. Firstly, we note that moving the 
function (ip) outside the entire electromagnetic-sector Lagrangian does not affect this issue. Indeed, multiplying 
the second term of (|21|) by any factor / affects both the fermionic kinetic term ipdfj^ip and the vertex ipAi^ip; however, 
the fermionic fleld enters quadratically in both expressions. After canonical normalization of the fermionic field, the 
factor / drops out from the physical value of the electric coupling constant. It is also difficult to imagine how one 
may try to modify the structure of the covariant derivative without spoiling gauge invariance. 

Secondly, one may try to arrange for a time evolution of / during reheating in such a way that Cphysicai is brought 
from a very tiny value at the end of inflation (so to avoid the strong coupling problem during inflation) to the present 
value before the onset of Big-Bang Nucleosynthesis (given that only a fractional discrepancy from the current value 
can be tolerated then [60|,|61j). We stress that this requires a huge change of /, which can be difficult to accomplish 
without disrupting the result for the magnetic field achieved during inflation. The comoving energy densities 

PB ^ TT^ I dkk^V'^= I dkpBk 



Pe = 

need to satisfy 



27r 

^ I dkP fv ~jv) = I dk -pEk (23) 



d _ _ r _ 

(pEk + PBk) ^ {PEk - PBk) (24) 
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If the electric component in this expression can be neglected, one finds pBk oc alternatively, if the magnetic 
component can be neglected, one finds pEk In general, achieving such a large change in / during reheating 

does not appear feasible. 

Thirdly, one may abandon the idea of identifying A^^ with the electromagnetic field, but still generate a large scale 
value of a hidden sector and weakly coupled A^, and then try to convert it to an electromagnetic field through some 
coupling. For instance, gauge invariance allows for 



P x-cm 



(25) 



The first coupling was originally proposed in [62l | , and one can promote x from a constant parameter to the expectation 
value of a scalar field; in the second case, x is a pseudo-scalar function. One could imagine that x experiences a quick 
transition from zero to a nonvanishing value x* at some given time r* after inflation, when / has set to one (we model 
the transition with a step function; clearly this approximation will break at very small scales) . Solving the equations 
of motion in vacuum at leading order in x* , and requiring continuity of the vector potentials at r* , one finds that the 
turning on of x results in a partial conversion of the "electric" or of the "magnetic" component of A^ into our electric 
field: f °™ ~ x* in the first case in (pS)) . and f °™ ~ x* S in the second case. This solution is obtained in absence of 
any charged particle. However, as soon as a plasma is formed, it shuts off the electric field well before it can convert 
into a magnetic field. If already present at , the plasma would prevent any electric field generation at all. 

In conclusion, none of these attempts appears to provide a solution to the strong coupling problem. We believe 
that a solution, if at all possible, will require a more radical modification of the model than those mentioned here. 

III. OVERVIEW OF THE MECHANISM AND PHENOMENOLOGY 

In this Section, we describe how the model impacts the cosmological perturbations. The emphasis is in describing 
the key points, and in summarizing the results. The actual rigorous computations are performed in Sections IIVI and 
lYl 



Summarized Discussion of the Mechanism 



In the last Section, we showed that the homogeneous condensate (po{t) = {(p{t,x)) leads to the production of a 
scale invariant spectrum of primordial "magnetic" fields via the coupling l'^{ip)F'^. However, this same coupling 
also provides a channel for the produced gauge field fluctuations to feed back into the perturbations of the inflaton 
dip{t,x) = ip{t,x) — (po{t). Heuristically, we can see this effect by looking at the equation of motion for the inflaton 
perturbations. This has the general form 



Sf" + 2H5lp' - ASip 



- B"^ 



(26) 



where % = a' /a, Ai^ <C H^, A = didi, and • • • denotes gravitational interactions and also terms involving more 
derivatives of /(</?), all of which give subdominant corrections to the phenomenology. We see from ([^51) that the large 
scale gauge field fiuctuations can source infiaton perturbations. Schematically, the solution of (pS)) behaves as 



S(p = difivac +£^sourccd 
homogeneous particular 



(27) 



The homogeneous solution is just the usual quantum vacuum fluctuation amplified by the quasi de Sitter expansion, 
Sfvac ^ H/{2tt). On the other hand, we have an additional contribution which is sourced by the gauge field fluc- 
tuations. This new sourced contribution is uncorrelated with the vacuum fluctuations, hence its contribution in the 
n-point correlation functions will add incoherently with the standard results. Moreover, the source contribution to 
5if is highly nongaussian; it is bilinear in the (nearly) gaussian gauge field fiuctuations. The curvature fluctuation 
Q = —M^S(p therefore will also be characterized by a new contribution that is nongaussian and uncorrelated with the 
vacuum part. 

A completely analogous discussion applies also to the tensor perturbations. Since the produced gauge fleld fluctua- 
tions carry anisotropic stress/energy, they provide a new (essentially classical) source of gravitational waves which is 
uncorrelated with the standard vacuum fluctuations. 

The underlying physics discussed here is very similar to the inverse decay processes that have been computed in [TtI - 
[igj and also the rescattering effects considered in [20l - |22| . However, in both of those cases the second order "sourcing" 
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of inflaton perturbations occurred near horizon crossing, leading to a (nearly) equilateral bispectrum. Here, on the 
other hand, we have a source term that is most significant on very large scales. This large scale entropy mode leads to 
a (logarithmic) time evolution of Csourccd on super-horizon scales, and consequently the bispectrum is of nearly local 
shape. 

B. Summary of the Key Phenomenology 

We define the power spectrum of curvature fluctuations via 

i(kCp) = ^Pdk)S^'^{k + k') (28) 

The final result for , evaluated on large scales and at the end of inflation, is 

Pc=V[l + 1927'7VgMB (^^tot - A^cmb)] (29) 

where P^/^ = ^^^^| is the amplitude of the power spectrum of the vacuum modes, iV^MB denotes the number of 
e-folds between the moment at which the large scales CMB modes leave the horizon and the end of inflation, and 
A^tot is the total number of e-folds of inflation. The flrst term in ([29)) is the usual contribution from the vacuum while 
the second term is the "sourced" contribution described heuristically in the last Subsection. Here we work to leading 
order in slow roll parameters, so the spectrum (j29p is exactly flat. In a more complete computation we would see 
small departures from scale invariance oc k"'~^. 

We can physically understand the structure of as follows. The suppression <C 1 in the second term arises 
simply because we are computing an effect which is higher order in perturbation theory.'^ The factors of A^cmb, on 
the other hand, arise due to the logarithmic time evolution, Csourccd ^ In a, outside the horizon. Such growth is 
consistent since we have large scale entropy perturbations playing an important role in the dynamics. Finally, the 
factor iVtot — NcMB is related to the phase space of contributing gauge field fluctuations. It is related to the number 
of B-modes that source the inflaton perturbation, and it is the counterpart of the logarithmic enhancement in the 
background density We explain this in details after eq. (I57|) . We should stress that equation ([29]) is valid only 

when TVtot ^ ^cmb- Otherwise the factor A^tot — Ncmb is replaced by an order one factor. 

Throughout this paper we will require that the sourced contribution to the power spectrum is subdominant: 

192 V N^MB {Ntot - Ncmb) < 1 (30) 

This yields a constraint on the total number of e-foldings iVtot — Nqmb < O{10~^)'P~^ (taking iVcMB 60 for 
illustration) which is considerably more stringent than the backreaction bounds discussed in the last Section. 
The bispectrum is given by the 3-point correlation function 

(Cfc, Cfc, CkJ ^ 5^'^ {ki + k2 + h) (31) 

We have found that our bispectrum is very close to the local shape. Indeed, the "cosine" between our bispectrum 
and the local shape, as defined in (63j (this is a measure on how well a template reproduced a given bi-spectrum) , is 
about 0.98. Nevertheless, here we will retain the full momentum dependence of the bispectrum, since it has a simple 
analytical shape. 

As is conventional in the literature (see e.g. [63]), we defined a /c-dependcnt nonlinearity parameter from computing 
the bispectrum and the power spectrum of C, and by comparing them with those obtained from 

c (f ) = + 1 /nl [c| {£) - ie, m] (32) 



^ See [27| for more discussion on the counting of such factors in models with particle production. 
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where Q is gaussian. An explicit computation gives the result 

l + COS^ (^1 , ^2 ) 1+COS^ (^1 , ^3 ) l-j-COS^ (^2 , ^^3^ 



r n \ /■equiv. local ^, 



3 kf k'^ ' k^ k^ ^ k^ k^ 



4 1,1 



fef"fc|' k'{ k'^ j,3 /j.3 
^cquiv. local ^ mOWaMBlA^tot-A^CMB) 

- (^^) (^tot-A^CMB) (33) 

where we have assumed that (l30l) is satisfied. We can a-posteriori see that this condition is indeed always satisfied 
whenever the result ([33]) is within the observational limits. 

We have defined our "equivalent" local nonlinearity parameter as follows: we first note that both our bispectrum 
and the local template are enhanced in the squeezed limit. In this limit, our bispectrum satisfies 



^1 (Cfe, Cfe, Cfc3> « 1 + C0S2 (^1, 4) , ks « fcl 



2 



(X COS e Yi? + sin eY^ , e = tan"^ ~ 0.22 (34) 
" 2V5 

where Y^'^ are normalized spherical harmonics, characterized by the angle between fci and k^. The average of /nl 
over all possible values of this angle is then equal to the av erag e of /^^'^' Note that the current CMB limit on 
local shape nongaussianity is -10 < fj^l""^ < 74 at 95% CL [ea]. 

Using the local template to study this signature is clearly a good approximation, given that the local template 
is characterized by the monopole part only in (p4l) . and that cos e ~ 0.98. However, we note more than 1/5 of 
the amplitude in p4p is contributed by the quadrupole part. This is in contrast to what typically happens in non- 
gaussianity from scalar fields only, where the quadrupole and higher harmonic terms in the squeezed limit expansion 
are suppressed by powers of k^/ki and give a negligible contribution in this limit [26|. The angular dependence is 
imprinted by a "directionality" generated by the largest wavelength mode fc^^, seen by the smaller modes when they 
leave the horizon. In the scalar case, the directionality is typically due to a gradient, and therefore it vanishes in the 
fcs — >■ limit. In our case, the directionality is due to the polarization of the vectors, and it therefore remains finite 
in the limit. 

Let us also discuss the parametrical dependence of /^^'^' The factor iV^j^g is due to the super-horizon growth 
of the three modes used in computing the three point function. The factor iVtot — .^cmb is due to the number of 
super- horizon modes that contribute to the correlator, analogously to what we have described in relation to ((29|) . Also 
in this case, this factor is replaced by an order one factor if iVtot is close to A^'cmb- 

The "magnetic" fields also produce gravity wave modes, which add incoherently with the vacuum ones. The power 
of gravity waves produced during inflation is conventionally parametrized by the ratio r of their power divided by the 
scalar power. We find 

r = ^ 16 c ^ + ' ^ ^CMB (^tot - Ncmb) , . 

1 + 192 V N^M^iNtot- Ncmb) ^ ' 

which gives the standard result r w 16e when the vacuum modes dominate the power spectrum of curvature fiuctua- 
tions. 



IV. SCALAR PERTURBATIONS 



To encode the effect of the gauge field on the cosmological perturbations we need to study the perturbations up to 
second order. Therefore, we decompose 

ip = LPQ + Sup + 52<p , gtj.v = gtj.va + Sig^ii, ^ 52gtiu , A^, (36) 

The gauge field has no zero order part, and, as we discuss below, we do not need to evaluate it at second order. 

Since the gauge field has no zero order part, the metric/inflaton perturbations do not mix with the gauge field 
modes at linear order; this is because the gauge field enters already quadratically in the action for the perturbations, 
through the expansion of the last term in This is the same reason that in the previous Section allowed us to 
compute disregarding infiaton and metric perturbations. We note that the gauge field can still affect the first 
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order metric/inflaton perturbations through its backreaction on the background evolution. This can be disregarded 
under the assumption that the two conditions (IT5|) and hold. 

Therefore, at the linearized level the standard results of scalar field inflation hold. We work in the spatially flat gauge 
for the scalar perturbations, Sgij = 0. In this gauge, the curvature perturbation on uniform density hypersurfaces is 

= —M- Sip. As we show in Appendix lAl one finds 



OT 



2 \ 1 



Slip = , 



(37) 



which is the standard equation for the Mukhanov-Sasaki [Gg, |63| variable rewritten in terms of C; the last term in this 
equation has been simplified using slow roll approximation. Finally, T~L = \ = o,H . Here we have disregarded correc- 
tions to the effective mass due to backreaction effects; see Section |lll In the end we will only perform computations 
at leading order in slow roll parameters, hence this neglect has no impact on our final results. 

To compute the effect of the gauge fields on the cosmological perturbations we expand all the equations of the 
system at second order in the perturbations. We combine these equations in the same formal way that they are 
combined to obtain ([37]). In this way we obtain a "master equation" for 52'p that does not contain any 52g^v mode. 
As we show in Appendix \^ this equation reads 



d 

2H—-A 

OT 



- 3 



Ml! 



S2tp = Ji [A^] + J2 (Siff , (Sigf , Slip X Sig 



(38) 



where 



2HM^ 



dr (a^V -{Ex B) 



(39) 



As for the linear theory equation, we have disregarded a backreaction-induced correction to the effective mass. 

Note that our model ^ contains higher dimension interactions between the infiaton and gauge fields of the form 
'^,i'^niSip)"'F'^ which arise from expanding the coupling function l'^{p) in powers of Sip. These couplings will enter 
into the calculation explicitly at higher order in perturbation theory. Using the in-in formalism, we have verified that 
such high dimension operators do not modify our leading order results for the spectrum and bispectrum. 

The right hand side of psp comprises of two sources for S2P] the first source contains terms at second order in the 
gauge perturbations, and has been completely given in p9l) . The second source contains terms that are the product 
of two first order infiaton perturbations, or of two first order metric perturbations (not only the scalar ones), or of 
one first order infiaton perturbations times one first order metric perturbation. We note that no "mix source" of the 
type Slip X A or of the type Sig^j^i, x in present in ([55)) . because does not enter linearly in (j?]). 

Expression was first obtained in [i^ by extremizing the cubic order action of the perturbations. This 

is equivalent to working directly with the equations expanded at second order, and we have verified that our result 



coincides with that of [43[ . The source J2 is the standard result obtained at second order in single scalar field inflation. 



The scalar part of this expression in the gauge we have adopted is explicitly given in [6S 
Therefore, the infiaton perturbation is formally given by 

Sip = {Slip + (52¥'|sourccd by Ja ) + '52<y5|sourccd by Ji 



(40) 



The part in parenthesis is the standard result obtained in single scalar field infiation, with only negligible corrections 
coming from the backreaction of the gauge field on the background dynamics. The 52'p|sourcod by Jj term is clearly 
negligible in the primordial power spectrum, and also leads to unobservable non-gaussianity |2|-|5|. We therefore 
disregard it in this work. The last term in (|40)) encodes the effect of the gauge field on the infiaton perturbations. 
We note that this term is uncorrelated with the other two, since the quantum/statistical operators entering in Ji are 
those of the gauge field. Therefore, we are interested in computing 



(Sip'^) ~ (| ((5 (/'vacuum)^ ^ + (^{Sipsomcedf'j , (^V^) - ^ (<5<Psourced)^ ^ 



(41) 



where 



Sipv 



Slip , Sipsouvccd = S2'p\ 



sourccd by Ji 



(42) 



We combine p7p and p8p in a unique equation for Sip = (5 (/'vacuum + '^V'sourced, where the two quantities are defined 
in (1421) . We approximate this equation in slow roll approximation and we keep only the leading source term that 
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arises from the direct I F coupling. This gives 



+ 2-dr - A 
a 



5lp ^ J 



^21 ^_ /I^ lE^-B^ 
J V e Af L 



(43) 



We note that the source — interacts with the inflaton perturbation with a strength that is gravitationahy 
suppressed but slow roll 1 / ^Je enhanced (this is one of the enhancements that make non-gaussianity visible in the 
model). The remaining terms in p9|) have the same scale dependence, but an interaction strength that is both 
gravitationahy and slow roll -^/e suppressed. The same suppression characterizes all the terms in J2. Therefore, the 
dominant source is that one arising from the direct inflaton-gauge field coupling P {(p)F'^. We explicitly see that 
simply computing the effect on S(p from the gauge fields in an unperturbed metric reproduces the leading results for 

In our calculation we estimate the curvature perturbation in flat gauge as C = —^6(p. In principle, there could 

be contributions to C at second order coming directly from the produced gauge field fluctuations (as opposed to 
the contribution from feed-back on the scalar inflaton) which could be estimated parametrically as ^ -^-^Spgaugc- 
During inflation such contributions will give corrections to the spectrum and bispectrum that are suppressed by large 
logarithms, with respect to what we compute here. Moreover, such "direct" contributions to C will be diluted during 
the reheating phase, when the equation of state of the dominating plasma is either the one of matter (for perturbative 
reheating), or intermediate between the one of matter and radiation [69^ (in the nonperturbative case). 

A. The two point correlation function, and the correction to the power spectrum 



We are interested in the primordial curvature perturbation, given by 

C{t,x) = -—S(p{t,x) 



(44) 



The two point correlation function in momentum space is related to the power spectrum by the standard expression 



a a 

where all quantities are evaluated at some time r, and where we use the convention 

1 r d^k 



Stp{T,x) 



a J (27r)3/2 "^fc 
We define the Fourier transform of the source as 



6ipj:{T) e**=-^ , Sifij: = S(t>k a (k^ + 5cj)l a) (-fc) 



-ik-x J 



(27r)^/^ 

where ([T2|) and (|43|) have been used. We find 

(r) f d^p 



d^p 



2eM, 



■E 



P A, A' 



4p J (27r) 



exip) ■ (k - p 



3/2 



En • Et ^ — En ■ Br -' 

y k—p y k—p 



(27r)^/^ 
ax (p) + a{ {~p) 



9AA' 



|pll/2|fc_pll/2 |p^3/2|;,_p^3/2 



a\' [k—p] + a 



(45) 



(46) 



(47) 



(48) 



We stress that the first term in the square parenthesis in the first line of (|48p is the ^^-contribution to the source, 
while the second term is the B— contribution. 

Combining eqs. P5|) . and (|T7|) . the equation for the Fourier modes of the inflaton perturbations reads 



(49) 



where all quantities are evaluated at the time t. 
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The homogeneous solution to this equation is the standard vacuum solution, which leads to the standard result 
(fro)) for the power spectrum (|^5|) . The homogeneous equation is obtained with the Green function method 



sourccd 



dr'Gfc (r, r') Jr:(T') 



(50) 



The leading order result in slow roll approximation is obtained by using the retarded Green function of eq. (|49|) in 
de Sitter space. 



Gfe (r, r') 



Using the source (|48)) . and the identity 



1 



3r 



- [kr' cos (fcr') - sin {kr')] , |fc t| < 1 
|fcT|, Ifcr'l < 1 



|eA (p) ■ ex' iq)\ 



(51) 



(52) 



we obtain 



a a 



I sourced 



H'5(^)[k + k') r r ^y' f dy" 



9eAf2 k^ 



1 + cos^ 



y' 



81 



|91 



k — q 



k — q 



18 cos 



k~ qj 



k — q 



(53) 



where we have introduced the dimensionless integration variables q = p/k, y' — —kr', and y" = —kr". The momentum 
integral need to be restricted so that the gauge modes participating in the original convolution were inside the horizon 
at the start of inflation (otherwise they would not be produced by this mechanism). This means 



k1 Ak-q\> 



1 



k iTin 



(54) 



The time integrations are instead restricted to times which are between Tin and r, and for which the sourcing modes 

have exited the horizon. This means Tin, — lir, — t^^''^ — < t' , t" < t. We do not need to include Tn in this condition 

' ipr \k-p\ ' 

thanks to (IMl) . Therefore, the time integrals in ([55)) are restricted to 



k\T\ < y', y" < Min 



1 



kT \k-q] 



(55) 



In Section |TI| we saw that the energy density in the produced field is logarithmically enhanced with respect to 
that in the field. Eq. ()53l) shows that the an analogous logarithmic enhancement takes place in the source of dip. 
In this case, all the three integrals in (|53p present an enhancement. Disregarding the subdominant .E— contribution, 
and performing the time integrals, the expression (|53p gives 



a a 



sourccd 



9ir_ 
7m 



(5(^) [k + k'^j r <fiq 1 + cos^ k-qj 



k^ 



Min 



In" 



\k-q\ 



k ~ q 



k\T\ 



(56) 



The momentum integral has most of its support at the two logarithmic poles; due to the symmetry between the two 
poles, we can simply evaluate the integral for \q\ ~ ^^—^ ^ 1 and multiply the result by two: 



12 H* S^^Uk + k') 1 

^7mI P ''''W\'"^'^"" 



(57) 
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The first logarithmic enhancement is due to the growth of the two modes and in the super-horizon regime. 
The growth is due to the presence of the entropy modes A^. At the end of inflation In^ = ^cmb- ^he second 
enhancement In (A:|rin|) — iVtot — -^cmb is due to the number of gauge field modes that source the infiaton perturbation. 
More specifically, we have seen that each sourced infiaton mode is obtained as a convolution of two gauge field modes. 
The momenta of these two modes need to add up to the momentum of the infiaton mode. The second enhancement 
occurs in the IR limit of one of the two gauge modes. This enhancement is the counterpart of the enhancement taking 
place for pB, which is also due to the number of large wavelength modes produced during inflation. 

As we discussed, the contribution ([57]) adds up incoherently with the vacuum one in ()45p . Using the relation 
C = ^J^^V, the sum gives 

Pc lend inflation ~ 7^ [l + 192 N^^B (^tot ^ ^VcMb)] (58) 

where we recall that V is the contribution from the vacuum term, for which the standard slow roll expression (jl9p 
holds. We also remind that iVcMB — 50 — 60 is the number of e-folds before the end of inflation when the largest 
scale CMB modes left the horizon, while iVtot is the total number of e-folds of inflation. The enhancement from the 
momentum integral takes place for iVtot ^ .^cmb; if inflation only lasted about the observed number of e-folds, then 
we estimate that the final momentum integral produces an order one result, so that the result (l58l) remains valid as 
an order of magnitude estimate. 

For TVtot ^ -^CMB ^ 60, the ratio between the sourced and the standard power spectrum is ~ 1.7 • 10^^ A'tot- The 
standard term dominates provided that infiation lasted less than about 600 e-folds. In the work, we assume that this 
is the case. 



B. The three point correlation function, and observable non-gaussianity 

We are interested in the three point correlation function of C as a measure of non-gaussianity. A common 
parametrization of nongaussianity is the nonlinearity parameter /atl, introduced by assuming that the curvature 
perturbation may be expanded as 

C (x) = (x) + I /nl [C, (x) ~ m (59) 

where Cgi^) is a gaussian random field (see our discussion in 18] for a detailed explanation of the sign conventions 
and 2tt factors in the following expressions). Both C and Q may be decomposed as in so that the relation between 
the q-modes of the Fourier decomposition is 

^fc = ^S.fe + ^ ^NL J ^^^Cg.kCg,k^p (60) 

The three point correlator of Cg vanishes, as this field is gaussian. However, due to the quadratic term in (I59p . the 
three point correlator of C is nonvanishing, and can be expressed through a sum of two point correlators of Cg- One 
finds 

(C£. Cs, (£3) - ^ {27rf' /nlPc ikf S^'^ {h + k, + k,) ^ (61) 

where the power spectrum was defined in (j45l) . To obtain this expression, one identifies the two point function of Q 
with that of C,g (as the difference is subleading in a perturbative expansion), and disregards the mild scale dependence 
of the power spectrum. 

By evaluating (C'^), and by inserting it in (j6ip . one defines an "effective" (momentum dependent) nonlinearity 
parameter, even when the intrinsic nongaussianity is not of the local form (j59l) . The dependence of /nl on the 
relative size of the momenta is denoted as "shape" of the non-gaussianity. Ref 63] provides a method to evaluate 
whether the shape obtained in a given model is well reproduced by the local template (namely, /nl constant in (|6ip) 
or by any other template employed in data analysis. 

As nongaussianity from the vacuum term is negligible, we need to compute 

r r r \~ /^^^^- ^^'^^ ((^9) 



where each expression is evaluated at some given time r. 
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We proceed as in the previous Subsection by inserting the source (|48)) into ([50]) and by evaluating the correlator. 
Keeping only the dominant "magnetic" source, we obtain 



n 



1 729 



irccd 



(r) 87r9/2 g3/2^j3 11 

XfJaiAaf^^^^ (fci -P1+P2) ^CT2A 



\p^?'^\h-p^\^/^ 
,5(3) (a?2 -P2+ Ps) S^.X.S^^'' (fcs - P3 +Pl 



where the integration regions are bounded in an analogous way as for the two point function: 



\Pt\, \k^-p^\ > -r-l ' I'^l < 1^*1 < Min 

Tin 



1 



1 



\p^\' \h-p,\ 

Performing the time integrals and employing the 5— fonctions, we obtain 
27 



(63) 



(64) 



87,9/2 £3/2^/3 



X / d^p In Min 



(fci + fc2 + ^3) 
1 1 



\p\\r\' |fci-:?^|r| 



In Min 



1 



1 



\ki~p\\T\ \k3+p\\T\ 



In Min 



\h+p\\r\' \P\M 



E., 4'^^* iP) 4''^ iP) E.. 4'-'* {P- ki) ef^' {p- fci) E.3 ^f'^* {P + k.) ei'^^ (p + fca) 



|pl3|p-fcl|3 \p + k3\ 



(65) 



where, due to (|64p . the integration region is delimited by 



IpI, |p-fci|, W+ksl > 



(66) 



The momentum integral in (j65p has most of its support at the three logarithmic poles; each pole occurs when one 
of the three quantities in (|66)) reaches its minimal value . Formally, 



n 



ki 



sourccd 



fci, k2, fca 



fcl, fc2, ks 



(67) 



where C refers to the contribution to the integral in (|55t from the region close to the pole indicated by the suffix. 

To evaluate the contribution from the second region, we redefine the integration variable as p — >■ p + fci; we then 
see that this contribution is formally equal to the contribution from the first region, provided the external momenta 
in the first region are changed as fci — >■ ^2, ^2 — J' ^3, and fca -> fci. Analogously, to evaluate the contribution from the 
third region, we redefine the integration variable as p — >■ p — fca; we then see that this contribution is formally equal 
to the contribution from the first region, provided the external momenta in the first region are changed as ki — >■ k^, 
^2 — > fcl, and fca — >■ ^2- We therefore have 



n 



fcl, fc2, fca 



permutations 



To evaluate this contribution, we use the identity 

4^^* = '^^J ^ Pi Pi 



(68) 



(69) 
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and we obtain 




sourccd 



27 

87,9/2 £3/2^/3 



6^^'^ (fcl + k2 + fcs) 



X < In , ^. . In , , In Min 
ki r fcg T 



1 + cos ( fcl , h 



1,31,3 



In Min [fci |ri„|, fca |rir 



+permutations 



(70) 



Assuming that the external momenta are not too hierarchical, we disregard the difference among them in the 
argument of the logarithms. The expression for the correlator, evaluated at the end of inflation, then simplifles to 



n 



ki 



sourccd 



144 \I^V^'^ (5(3) (fcl + fc2 + 4) iV^j^B (^tot - iVcMB) 



1 + cos' 



(^1, ^2) 



1 + COS' 



(fcl, ^3) 



1 + COS' 



3 i,3 
2 



^3 



1,3 1,3 
'^1 '''3 



+ 



(fc2, 



)l3 l3 
"■2 "-3 



(71) 



We can now evaluate the 3-point function of the curvature perturbation using the relation C = —-^Slp and introduce 
a momentum-dependent nonlinearity parameter by comparison with (I6ip . We find: 



l+cos^(fci, /C2) 1+003^(^1,^3) lH-cos^(fc2, /i^a) 



■cquiv. local 



/•cquiv. local 
JNL 



1280 J^N^^s (iVtot - A^cmb) 



(72) 



Pc(fc)2- 

If the power spectrum is dominated by the vacuum fluctuations (which we assume for this work) then we have 

j^cqui 
INL 



Pc(fc) « V and rZ"- ~ USOVN^M^ (Mot - ^Vcmb). 



V. TENSOR MODES 

Production of gauge field fluctuations during inflation and its effect on the curvature perturbation have been 
discussed in the previous Sections. The produced gauge quanta, however, affect not only the scalar but also tensor 
perturbations (gravity waves) . Metric perturbations couple to each content of the universe and are inevitably sourced 
by the produced gauge fleld fluctuations. To see the effect, we consider the transverse and traceless components of 



the spatial metric perturbations: 



a {5ij + hij ) , with /i^; 



and dihij 



0. As the matter content is scalar and 



vector, hij is the only tensor perturbations in the model. 

From the Einstein equations, one finds the same equations for hij in terms of the physical Ei and Bi as in [IE 



1 

2^ 



2^9. 
a 



A h, 



— [E^Ej + BiBj)^'^ 



(73) 



where TT denotes the transverse and traceless projection of the spatial components of the energy-momentum tensor 
of the gauge field. ^ Since the gauge field has no expectation value, there is no coupling to the tensor perturbations 
at linearized level. Thus ([75)1 is in fact up to second order, and the right-hand side should in principle contain the 
source terms from squares of the first-order inflaton and metric perturbations. As for 82'^, such source terms are 
uncorrelated and subdominant to those from the gauge field. 



This projection can be done by an operator Ojj 
\Pij{k)Pim(k), where Pij(k) = 5ij - kikj . 



which can be expressed in the momentum space as Oij i^(k) = Pii{k)Pjm\k) — 
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Tensor modes (or GW) are transverse and traceless part of Sgij and have two physical degrees of freedom, the 
left-handed (L) and right-handed (i?). It is convenient to decompose tensor perturbations as 



hij (t, x) = 



(27r 



,3/2 



s=L,R 



ny-,(fc)/i,(fc) 



(74) 



"''^^ (k) ef (k) , which 



where hs{k) = hs(k)as{k) + /i* (fc)aj(— A:), and the helicity projectors are Ilij,L/R{k) = e^^'' (fc) ej^' ' 

clearly have the properties Hu^s (k) = h Ilij,s {k) = 0. Note that hg (k) depends only on the magnitude of k. Since the 
mechanism of GW production is analogous to that of the curvature perturbations, presented in detail in the previous 
Sections, we merely show the result of our computation here. Define GW power spectrum Pl/r. hi the usual way. 



Psik)'-^S^'\k- 



k'] 



(75) 



where s = L/R. As for the scalar perturbations, the GW modes produced by the gauge quanta are uncorrelated with 
those from vacuum fluctuations, and so the two contributions simply add up in the power spectrum. The two helicity 
states are produced in the same amount, and their sum gives 



Pgw (k) 



2IP_ 

7r2M2 



1 



6IP_ 



In" 



In- 



H 



(76) 



where the first term in the square brackets is the contribution from the vacuum and the latter from the source. 
Evaluating this expression at the end of infiation gives 



iiflation 



2H^ 



1 



N, 



CMB 



CMB 



(77) 



As expected, the standard vacuum part dominates in the regime we are interested in: when the vacuum contribution 
to is dominant over that from the source, the same is also true for the GW power spectrum; this is due to the 
fact that the tensor modes are only produced gravitationally, while the dominant source of the scalar modes is the 
direct inflaton-gauge field coupling; (this coupling is mathematically enhanced with respect to the gravitational one 
by 1/e). In this regime, the tensor-to-scalar ratio r reproduces the standard result:^ 



16e 



(78) 



VI. VALIDITY OF PERTURBATION THEORY 



It is important to verify that cosmological fluctuations remain small in our model, to justify a perturbative analysis. 
In this Section we consider several different diagnostics, arguing that perturbation theory is well under control, and 
that the second order results we have computed are subdominant with respect to the standard flrst order results 
(apart of course for what concerns the non-gaussian nature of C, which is not present in the linearized theory). 

As discussed previously, the produced gauge fleld fluctuations backreact on the classical background, by contributing 
to the total energy density in the Friedmann equation and by introducing dissipation into the equation for the 
homogeneous inflaton. A first consistency check is to ensure that these effects are negligible, which we have already 
shown in Section HIl 

Another important diagnostic is the variance of curvature fluctuations, (C^)^^'^ (in general, we need to compute the 
variance of a perturbation as a measure of the value of that perturbation in real space) . Using our previous result for 
the 2-point correlation function in momentum space, it is straightforward to compute the total variance of curvature 
fluctuations 



(C^(r,a-)) = {e{r,x))\+ (C'(r,f))| 



sourccd 



V 



dk 
~k 

a{r) 



- 192V 
1 + 16Pln 



dk 



■In' 



k 

3 I «(t) 
ain 



a{T)H 



In 



a-niH 



(79) 



^ We define the tensor-to-scalar ratio in the usual way, by normalizing the power in GW to that in curvature perturbations: r = Pgw/Pc = 
{Pl+Pr)/Pc- 
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To compute (|79p we have integrated from fcmin = ainff to fcmax = o,{t)H , since we should count the phase space of 
super-horizon modes from the onset of inflation to the time r. Equation ([7^ shows that (C^)^^^ ^ 1 generically and, 
moreover, contributions from sourcing effects are controlled by the small parameter V ^ 10^^. 

Since the spectrum of sourced curvature fluctuations exhibits a logarithmic growth, one might worry about es- 
timators of inhomogeneity that contain derivatives of the fluctuations, for example the energy density or the Weyl 
curvature. We consider such diagnostic in the next two Subsections, which are devoted to the contributions from 
tensor and scalar modes, respectively. In the final Subsection we conclude that perturbation theory is well under 
control. 



A. Tensor Modes 



We first consider the contribution to the energy density and Weyl tensor from tensor perturbations, since they are 
technically simpler to evaluate. 

The energy density in gravitational waves contains a contribution due to the usual vacuum fiuctuations, and an 
uncorrelated contribution from sourcing effects. The former is given by 



5p~ 



GWl 



167r2 

t2 



The sourced contribution is given by 

(^PGW I sourced ~ 



47r2 a2 



dk 

T 



6g2 



In 



In' 



a{r) 



k^ In^ 



(80) 



a{T)H 



(81) 



In both ([80l) and ([8T|) we have integrated from kmin = ainH to fcmax = (^{t)H, as discussed above. We see that the 
sourced contribution is subdominant to the vacuum one, being suppressed by the small number /M^ ^ 10^^°. The 
constraint (Jpcw ^ SH'^Mp is easily satisfied. 

It is interesting to compare the time dependence in the final results ((80|) and (|8T|) with what would be obtained by 
inspecting the spectral densities. In position space, the sourced contribution exhibits a logarithmic growth, relative 
to the vacuum contribution (we note that, however, it remains always very subdominant). For the spectral densities, 
on the other hand, we see that the results differ by a power-law dependence on t (this does not imply that either 
result is growing with time; indeed, neither term grows, due to the time dependence of the scale factor outside the 
integral). This difference in behavior is not seen in the integrated result, when the phase space of relevant modes is 
accounted for. 

The Weyl tensor is a gauge invariant object which vanishes for an FRW universe and thus may be considered as a 
measure of inhomogeneities. A good estimator for the size of typical entries in the Weyl tensor (from gravity wave 
perturbations) is the quantity 



C^(r,f) = h" + Ah. 



(82) 



Refs. j7(j] used as a dimensionless diagnostic of the validity of perturbation theory the ratio between the entries of 
this tensor and those of the Ricci tensor, whose background components are of the order H^. Proceeding as above, 
the contribution from vacuum fiuctuations is 





r dk 


'm^k^T^' 


1 vac 


2^ J T 





For the sourced contribution we instead have 



sourced 
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dk 
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(83) 



(84) 
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Again, we see that the contribution from sourcing effects exhibits a logarithmic growth, but remains safely sub- 
dominant to the vacuum part. As for the energy density, the power-law time dependence in the spectral density is 
compensated by phase space factors. 



B. Scalar Modes 



We now consider contributions to the energy density and Weyl tensor from scalar fluctuations. The energy in 
inflaton fluctuations, Sp^p = ^{{^"^'Y + {^^v)^ + given by 



5p, 
H 
8^ 



^ I vac 
4 



-f Sp, 

1 + ST^ln 



V'Isourccd ' 

a(r) 



air) 



?ln^ 



air) 



(85) 



where the first term is the usual contribution from the vacuum, the second term arises due to the feed-back of produced 
gauge fluctuations into the inflaton perturbations, and r] « V^^/{3H^) is the standard slow roll parameter. We see 
that the condition Sp^ <C 3H^Mp is easily satisfied, and that the sourced contribution is subdominant. 

A good estimator for the size of the contribution of the scalar fluctuations to the Weyl tensor is the quantity 



(86) 



where $ and ^ are the gauge invariant Bardeen potentials. This estimator is discussed in some detail in Appendix 
A. Again, we compare this to TL^ to have a dimensionless diagnostic of the validity of perturbation theory. The 
contribution to our estimator (|86p from the usual vacuum fluctuations arises at linear order in perturbation theory 



C^(T,f)| 



^cr-ACi) 



(87) 



It is straightforward to compute the variance 



(88) 



where in the second line we have introduced the notation ^^(t) = —iir^^'^'P^^'^ (—t)^^''^ H!^^^]^{—kT) for the mode func- 
tions. Before moving on, it is worth commenting on the constancy of the result (j88p . Even if the Fourier transform 
of /SC,/'H? is decreasing in time in the super-horizon regime, the momentum integral that gives the variance (1881) is 
dominated by the UV cut-off k = aH] this phase space factor compensates for the factor T-L^^ — {aH)~^ . (Note that 
this same phase space compensation was seen also in Section |lT] when we computed the energy density in electric and 
magnetic flelds, showing that pE is "only" logarithmically suppressed with respect to pB^ even though the electric 
held Fourier modes decay as outside the horizon.) Despite not decreasing in time, the integrated result is ^ 1. 

At second order in perturbation theory there arises also a contribution to the Weyl tensor due to the feed-back of 
the produced gauge fluctuations on the inflaton and scalar metric perturbations. In Appendix A we show that this is 
of the form 



I sourced 



6(C^' - AC2) + T . 



(89) 



The source term T is quadratic in gauge fleld fluctuations and the full expression (given in Appendix A) is quite 
cumbersome. For our purposes, it sufflces to consider a representative contribution to the source term 



M2 



(90) 



During inflation, such quantities are tiny as compared to TH? by virtue of the backreaction constraint pe + Pb <SC 
3H^Mp. Moreover, notice that the source term T will quickly decay as after inflation has ended (at which point 



the production of gauge fluctuations comes to a halt and we recover the usual behavior pE ^ Pb 



we expect that ^2 becomes frozen on large scales after inflation, 
irrelevant for any late-time observable. 



Therefore the source term T in 



). By contrast, 
should be 
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Let us now compute the variance of the first term in (|89p . Following very closely our previous computations we find 



?£-M(C^(t,x))') ^96e'r^ln'{^]+--- (91) 

sourccd \ dii 



This exhibits a logarithmic growth during inflation but it always remains much smaller than the vacuum contribution. 



C. Summary 

To summarize, in this Section we have computed the contribution to various physical quantities from both the usual 
vacuum fluctuations and also from the feed-back of produced gauge field fluctuations on the inflaton and scalar metric 
perturbations. We have found that the "new" contributions from sourcing effects generically exhibit logarithmic 
growth during inflation and remain safely sub-dominant with respect to the standard result. This statement applies 
to the curvature perturbation, the energy density in fluctuations, and the Weyl tensor. We have computed our 
diagnostics during inflation. In principle, these results should be evolved through reheating and into the radiation 
phase, to ensure that the universe remains homogeneous and isotropic after the end of inflation. We do not expect 
any problem, since all the physical quantities at the end of inflation are dominated by the usual vacuum contribution 
(which does not lead to any unacceptable instability), since the energy density in the gauge field is subdominant, and 
since we assume a standard y/—gF^ term from the end of inflation onwards. 

Refs. [tO, [zH argued that the gauge field production in these frameworks leads to an unacceptable breakdown of 
homogeneity and isotropy soon after inflation, after the perturbations are matched to the post-inflationary epoch. We 
believe that the core of their claim is in the time dependent behavior of the spectral density of the sourced part of the 
Weyl tensor during inflation; they pointed out this feature in the scalar sector, but we performed computations also 
in the tensor sector because the identical feature appears also there, and the computations are simpler to present. We 
note that indeed there is a different time dependence in the spectral densities of eq. (|5(7|) vs. (|5T|) and of eq. (|83p vs. 
(|84p. We note, however, that this does not imply that the sourced term dominates over the vacuum one in physical 
integrated quantities like the energy density and the Weyl tensor (if the dominance does not take place before the 
end of inflation, it will not take place afterwards). All the sourced quantities that we have studied remain small both 
in real and momentum space, and do not grow as power law of the scale factor. 

VII. CONCLUSION 

In this paper, we have considered a simple model where the scalar inflaton is coupled to some U{1) gauge field in a 
way that would be typical for moduli or dilaton-like fields. We have seen that the time dependence of the condensate 
during infiation leads to a production of large scale gauge field fluctuations, analogous to the usual mechanism that 
amplifies the quantum vacuum fluctuations of the inflaton or gravity wave perturbations. Focusing on the case where 
the spectrum of produced "magnetic" fields is scale invariant, we have shown that (nearly) local nongaussianity is 
very naturally generated at the level /nl ^ 0(10), which should be probed in the near future. Although we have 
neglected slow roll corrections to the running of the spectrum and bispectrum, we can still see a logarithmic running 
of the effective Jnl parameter with scale. This arises since the "sourced" contribution to the curvature perturbation 
experiences a super-horizon evolution during inflation, due to the presence of large scale iso-curvature perturbations. 
Logarithmic running of this type may be of observational interest; see [72 74] for example. Moreover, since the 
nongaussian part of ^ is uncorrelated with the gaussian part, we have a non- hierarchical scaling which can lead to 
interesting signatures in probes that are sensitive to the global structure of the PDF . 

A novel feature of our result is the dependence of Jn l on the TVtot — Ncmb which measures the number of "extra" 
e-foldings of infiation, beyond the minimal A^cmb ^ 60 e-foldings between the end of inflation and horizon exit for 
CMB scales. Ordinarily, one would expect that such "extra" e-foldings are completely unobservable, since scalar 
modes which leave the horizon prior to A^cmb should just be absorbed into a renormalization of the homogeneous 
background. However, in our model the total duration of the quasi de Sitter phase impacts the energy density of 



produced gauge field [4j]. Indeed, if A'tot is too large then the backreaction of produced gauge fields will become 
appreciable and spoil inflation. In a restricted sense, we see that in this model the "extra" e-foldings of inflation can 
influence physical observables. 

We have seen that the correlation functions of scalar and tensor cosmological perturbations exhibit a logarithmic 
time dependence which is related to the growing phase space of produced gauge fluctuations. These logarithms should 
not be confused with the IR logs that have been discussed extensively in the literature in association with loop effects 
during inflation [75l - [78j (see also for a related discussion). In the case at hand the interpretation of this logarithmic 
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time dependence is straightforward. The production of gauge field fluctuations in our model arises simply because 
the effective gauge coupling is time dependent. This time dependence leads to a growth in the energy density of 
gauge fluctuations, which is drained from the scalar condensate. The energy transfer provides a physical clock, and 
its logarithmic growth is a real physical effect that is not related to the de Sitter background or to the quantization 
of gravitational fluctuations. 

There are a number of interesting directions for future research. It may be for example interesting to study effects 
that can arise if the vector field has a vacuum expectation value already at zeroth order in perturbation theory. It 
would also be interesting to consider different choices of coupling function, for example n 7^ 2, — 2. In connection with 
this, it would be interesting to explore in more detail the nongaussian phenomenology of our model, in particular 
as regards higher moments, scale dependence, and LSS probes. It may also be possible to obtain interesting signals 
for gravitational waves at interferometers. Finally, it would be worth to investigate whether this mechanism for the 
amp lification of the gauge "magnetic" modes can be consistently modified so to avoid the strong coupling problem of 
[44| . and therefore be used as a mechanism for magnetogenesis. 
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Appendix A: Cosmological Perturbation Theory 



We work in the spatially flat gauge and expand the metric up to second order in perturbation theory as 

300 = -a2(l + 2</.i +2<?!)2), 
go, = a^d,{Bi + B2) , 
gi-j = a^(5,y . 

Similarly, the scalar and gauge fields are expanded as 

^p{t,x) = (fioit) + Siip{t,x) + S2(p{t,x) , 
A^,{t,x) = {0,6iA,it,x)+d2A^{t,x)) . 

The equation of motion for the inflaton field is 

1 



(Al) 
(A2) 
(A3) 



(A4) 
(A5) 



(A6) 



We expand (jA6p up to second order, using the Einstein constraint equations to eliminate the metric fluctuations in 
order to close the system. This procedure has already been detailed in ; here we generalize those results to include 
also the presence of the gauge field. (The analogous procedure for an axial coupling is discussed in [isj.) 
At linear order in perturbations, the Klein-Gordon equation (jA6l) gives 

(A7) 



Slip" + 2'H5iip' - Adip + a^V^^^Siip — i^q^'i + 2a^l/y0i — ip^ABi = . 
From the — and Q — i Einstein equations we can derive the results 
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2M2 



[2a'^V4>i + ip'Q5iip' + a^V^^5i<p\ = 0, 



mi 



2Af2 
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Using these to eliminate 4>i and Bi from (jA7l) we obtain 



— 2-H — 



,(^0) 



M2\ 1 



= 0, 



(A8) 
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where we have made use of the background equations and work to leading order in slow roll parameters. 

At second order in perturbation theory the procedure is nearly identical, save for the appearance of source terms, 
which are quadratic in first order fluctuations. From the Klein- Gordon equation (jA6P we have 



//., 



[SiA'.SiA'^ - ^^6lAJ^^SlAJ + d^SiA,dj6iA^] 



(All) 
(A12) 



where • • • denotes terms of the form ((5i</?)^ which will not be important for our computation. These have already been 
studied in ^68:] and they are known to give a negligible contribution to nongaussianity The relevant constraint 
equations are 



nAB2 + ^ [2aV02 + ip'oS2ip' + aV^^2¥'] = S^^^ 



where the source terms are given explicitly by 

/2 



4a2M2 



[SiA'^SiA'^ + d^SiA,{d.,SiAj - djSiA.,)] + 



2a2M2 
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Again, we have suppressed terms of the form (Sicp)'^. Combining the constraints (jA13l) and (jA14[) with (lAlip we can 
eventually arrive at the master equation 
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^eijkdjSiAk respectively. Equation 



where the physical "electric" and "magnetic" fields are Ei = —^SiA'^ and B 
(UT8)) is the main result of this appendix. 

Before finish, we briefly consider the Weyl tensor, which is discussed in Section IVll The Weyl tensor from scalar 
perturbations has only an "electric" part which is given by 



(A19) 



where <i> and are gauge invariant potentials, defined explicitly in [79l|. A good estimator for the typical size of 
entries in (jAl9p is the quantity 



= A($ + ^) = A(P+ AB' , 



(A20) 



where in the second equality we have restricted to the flat slicing. 

Equation (IA20p can be re-written in a suggestive manner. Introducing the gauge invariant curvature perturbation 
(^n — — ^(S„(/5 the constraint equations can be written in a convenient form 



ABi = eCi , AB2 - eC2 + H-^S^^^ - (3 - e)^^^^) , 



(A21) 
(A22) 



at linear and quadratic order. It is now straightforward to expand (|A20p in perturbation theory. Working to leading 
order in slow roll parameters we find 



Cf = e(Cr-ACi), 

Cf ^ e(C2 - AC2) - (3 - 6)9,5(3) + n-^AS^^'> + n-^drS^^^ - (1 - e)5(2) 



(A23) 
(A24) 
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Appendix B: The In-In Formalism 



In this Appendix we verify that the Green's function method employed in the text is equivalent to the "in-in" 
formalism at leading order. To compute correlation functions using the in-in method, we must first identify the 
interaction Hamiltonian. To this end, we employ the Arnowitt-Deser-Misner (ADM) form of the metric and integrate 
out the lapse function and shift vector. This procedure has been described in [i^ and also [ill, so we do not reproduce 
the details of the calculation here. The quadratic action for the inflaton perturbation is 



S2 



^ / drd^xa^ 



' M2 



The cubic interaction terms in the Lagrangian are 



dTd^xa^^j-dLp 



- 



(Bl) 



(B2) 



Here we have suppressed terms of the form (Sip)^, which are irrelevant for our calculation. Varying (jBl|) and (jB2p 
and expanding Sip = Slip + S2(p reproduces exactly the master equation (jA18|) . 

For the choice of I{p) under consideration the interactions on the first line of (jB2l) are controlled by the dimensionful 



- jg-. In contrast, the interactions in the second line of (|B2[) are controlled by a coupling 



coupling 

Therefore in the slow roll limit, e ^ 1, the interactions on the first line of (|B2[) are the dominant ones. 
At leading order in a slow roll expansion, the interaction Hamiltonian Hj(t) ~ ~ J d^xa^Cs can be written as 



Hiit) 
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where C = —-^Sip is the curvature perturbation and the source Jk{T) was defined in (j47p . The in-in formula is 
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A key simplification arises from noting that the mode functions of the produced gauge fluctuations are real-valued, 
up to an irrelevant constant phase. This implies that the produced gauge field fluctuations are commuting variables, 
to a very good approximation. We have 



(B5) 



where it is understood that only super horizon modes are relevant. Consequently, the source terms Jq(t) in (jB4p are 
mutually commuting and they may be pulled out of the nested commutator. The remaining commutators are easily 
evaluated using the formula 



C£,(ri),C£,(r2) 



Gfc,(ri,T2) 
' ipl a{Ti)a{T2) 
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where the Green function was defined in ()5ip . This formula is valid only for ri > T2. 

Using the commutativity of the source terms and the formula (IB6|) it is straightforward to evaluate the sourced 
contribution to the n-point correlation functions of For the 2-point and 3-point we have 



(Cfe,C,-,(r)) 



sourced 



atpoj 



sourced 



dTidT2 Gfei (r, Ti)Gfe2 (t, T2) ( J^^ (^O^/fe^ (T2)) , 

> 

—] / dTldT2dTzGk^{T,Tl)Gk2{T,T2)Gk-AT,Tz) 



24 



These coincide exactly with what we obtained previously using the Green function method. We have also verified 
that the cross-correlation of gauge field fluctuations with the curvature perturbation agrees with what was presented 
in [iHl, at leading order. 
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